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1 Introduction 



This talk discusses recent developments in the theory of inclusive radiative B decays, 
where recent is defined as being after the CKM 2008 workshop. We focus mainly on 
non-perturbative aspects of i? — )■ since this is arguably the most important new 
development in the field. 

Inclusive radiative B decays, namely B — )■ X^'y, B — )■ ^^7, and B — t- X^l^l' 
are important probes of new physics. Since there are no tree-level fiavor changing 
neutral currents in the Standard Model (SM), these processes can only proceed via 
loop suppressed transitions. These loops are sensitive to the masses and coupling 
of heavier SM particles, but they can be also sensitive to masses and couplings of 
particles that appear in extensions of the SM. By having a good theoretical control 
of the SM prediction and in particular the non-perturbative effects, one can place 
important constraints on these models of "new physics" . Currently the experimental 
value of Bi{B X,7, > 1.6 GeV) = (3.55 ± 0.24 ± 0.09) ■ 10"^ P is in agreement 
with the three theoretical predictions, which differ in the way they address the photon 
cut effects. The first, Br(B ^ X,7, > 1.6 GeV) = (3.15±0.23) ■ 10"^ |2] by Misiak 
et al., assumes that at the value of E^ > 1.6 GeV one can ignore these cut effects. 
The second, Br(S ^ X,7, E^ > 1.6 GeV) = (2.98 ± 0.26) • 10"^ [3j by Becher and 
Neubert, is obtained by an MSOPE based analysis of the cut related effects. The 
third, Br(i? ^ X,7, E^ > l.Q GeV) = (3.47 ± 0.48 ± 0.17) • 10"^ |4| by Andersen and 
Gardi, is obtained by DGE based analysis of the cut related effects (We have added 
the second error of 5% from non-perturbative effects which was not included in |3]). 

Since the 6-quark mass, rrih, is much larger than Aqcd, one can expect that the 
partonic rate r(6 — )■ 57) is equal to r(_B — )■ ^^7), up to effects suppressed by some 
power of AqcB/^b- For a long time it was assumed that just like other inclusive 
B decays, such as B ^ XJu, these effects arise only at the second power of this 
ratio. We now know that it is not the case. In fact, non-perturbative effects arise 
already at order AQCD/''^fe- The reason is that unlike B — )• XJu, there is more than 
one operator in the effective Hamiltonian that can contribute to the decay. As it 
is well known, apart from the dipole operator (^77 = (— e/87r^)m5Scr^,,F^''(l + 75)6 
the operators Qsg = {—g/9>-K'^)mbS(7^yG^^{l + 75)6 and Q\ = {cb)v-A{sc)v-A are 
also importanci (see [5] for notation). Loosely speaking, we can convert the gluon 
or the quark-pair into a photon, but at a cost of a perturbative factor of or a 
non-perturbative factor of Aqcd/'"^6- While the former effects are well studied, the 
latter have received much less attention in the literature. It is the contribution of 
operators other than Q-jj that give rise to these enhanced non-perturbative effects. 
Furthermore, unlike non-perturbative contribution familiar from other inclusive B 
decays, these non-perturbative effects are not given in terms of matrix elements of 
local operators, but instead in terms of non-local ones. 

^Penguin operators contributions are suppressed by their small Wilson coefficients. 
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Taking these effects into account, the following picture emerges. In the endpoint 
region of rrih — 2E^ ~ Aqcd, the photon spectrum can be factorized, symbolically, as 

dT{B ^ X,-f) = H-J(g>S + H-J(g>s(^J+H-J(g>s(^J(^J. (1) 

The first term of ([1]) is the "direct photon" contribution familiar from the factor- 
ization formula for B — )■ Xul v in the endpoint region. The direct photon contribution 
arises from diagrams in which the photon couples directly to the weak vertex. H are 
hard functions parameterizing physics at the scale m^, J are jet functions describing 
the physics of the hadronic final state X with invariant mass Mx ~ ^J^n^K(^^^ and 
5* are soft functions incorporating hadronic physics associated with the scale Aqcd- 

The second and third terms corresponds to the "resolved" photon contribution in 
which the photon couples to light partons instead of coupling directly to the weak 
vertex. Unlike the direct photon contribution, the resolved photon contribution arises 
only at order Aqcd/'"^6 and higher. The second (third) term of ([1]) corresponds to 
the interference of an amplitude in which the photon does not couple directly to the 
weak vertex with an amplitude in which the photon couples (does not couple) directly 
to the weak vertex. The new jet functions J probe the hadronic substructure of the 
photon at a scale of order -yZEyVqcD- The new soft functions, s, describe the soft 
interactions between the hadronic substructure of the photon and other soft particles, 
e.g. the heavy quarks. Unlike the soft functions in the direct photon contribution, 
they contain non-localities in two light-cone directions. See [5] for a more detailed 
discussion. 

As one integrates over large enough portions of phase space, the direct photon con- 
tribution reduces to matrix elements of local operators multiplied by calculable short 
distance coefficients. The resolved photon contribution does not reduce to matrix 
elements of local operators, but instead to non-local matrix elements convoluted with 
calculable short distance coefficients. The resolved photon contribution arises only 
at order AgcD/'^b- Since non-local matrix elements are functions, they are harder 
to model. This ultimately limits the precision one can achieve in the theoretical pre- 
diction for inclusive radiative B decays and in using them as a tool to constrain new 
physics. 

The rest of the talk is organized as follows. In section |2] we discuss recent develop- 
ments that has to do with perturbative aspects of -B — > Xg'^. In section [3] we discuss 
recent developments that has to do with non-perturbative aspects of 5 — > X^'y. These 
have implications also for B — )■ X^pi which we comment on in section HI In section 
|5]we briefly review recent developments in 5 — )■ Xsl^l~. We present our conclusions 
in section [6l 

Two other issues that we do not have time to discuss in detail are resummation 
and extrapolation. The first has to do with photon energy cut related effects. We refer 
the reader to E. Gardi's talk at CKM 2008 [6] . The second is the fact that theory and 
experiment are not compared directly. The experimental value of |T] is extrapolated 
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from the measured ~ 1.9 GeV to > 1.6 GeV using [7J. Considering our 
new knowledge of non-perturbative effects and the newer Belle measurement with 
E^ > 1.7 GeV [Sj, which can potentially allow us to compare these extrapolation 
factors against data, it is time to revisit the issue. 

2 Recent developments in ^ ^ : 
Perturbative Aspects 

Recent developments in the understanding of of the perturbative aspects of -B — )■ Xg'y 
can be divided into three parts. The continued effort to complete the NNLO predic- 
tion for r(6 — i- S7), the calculation of the {AqcB/^b}^ corrections to the Qj^ — Qj^ 
contribution to the integrated rate at order a^, and the calculation of the subleading 
jet functions contribution to Qi^ — Q-jry contribution to the spectrum. In other words, 
these are all parts of the direct photon contribution to the integrated rate. We discuss 
each of these in turn. 

2.1 T{b S7) at NNLO 

As was discussed in the introduction, V{B — )■ = V{h — )■ 37) + O {Kqq-c / mij) . 

There is a continuing effort to complete the calculation of r(6 — )■ 57) at NNLO. For 
details we refer the reader to C. Greub talk at the CKM 2008 [9]. Out of the three 
necessary ingredients: matching at /i ~ M^j/, running from /i ~ My^/ to /i ~ m;,, and 
calculation of the matrix elements at 0{a1), the first two are complete and the third 
is almost done. 

Since CKM 2008 there has been several new calculations that aim to complete 
the last ingredient. The 0{a1) corrections from Qj^ — Qsg were calculated in [lOj . 
The authors conclude that this contribution "..will not alter the central value of [2] 
by more than 1%." Details on evaluation of the NNLO QCD corrections in the heavy 
charm limit {rric ^ rni,/2) appeared in [TT]. After the CKM 2010 workshop, the 
O{(3oa1) corrections from Qgg — Qgg were calculated in [12]. The authors conclude 
that the correction to the branching ratio "amounts to a relative shift of +0.12%" 
for the cut of E^ > 1.6 GeV. The results of [12] were confirmed in [13], which also 
calculated the O{f3oal) corrections from Qi — Qsg and Q2 — Qsg- The authors conclude 
that "numerical effects of all these quantities on the branching ratio remain within 
the ±3% perturbative uncertainty estimated in [2]"- Finally, as discussed in [llj . 
the complete 0{a1) calculation of Qi — Qi^ and Q2 — Qr--/ is underway, where the 
goal is to "make the perturbative uncertainties... negligible with respect to the non- 
perturbative... and experimental... ones". 
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2.2 C>(AQ(^j-,/m^) corrections to Tjy 



Considering only Q-j^ — Qj^, power corrections to the integrated rate are from order 
^QCD/^^b suppressed local operators. There are two such possible operators: the 
"kinetic" operator and the "chromomagnetic" operator. The perturbative coefficients 
of these operators were calculated at tree level long time ago in pH]. Recently, they 
were calculated at 0{as) in [15], where the authors presented analytical expressions 
for these coefficients. As for the numerical impact, they conclude that "The effect on 
the B X,7 rate is below 1% for < 1.8 GeV." 

2.3 Sub leading Jet Function contribution to dTjj 

Considering only the Q^^ — Q^^ contribution to the spectrum, one has the following 
symbolic factorization formula in the endpoint region: 

1 1 /A^ \ 

dTrr r^H-J(g)S+ — YH-J(g)Si + — yH-ji^S+Oi (2) 
mb ^-^ mb ^-^ \ J 

The factorization of the leading power term, H ■ J ® S, was proven in fl6\ [T7] . The 
factorization of the one of the possible subleading power terms J2i H ■ J Si, namely 
the subleading shape functions (SSF), was proven in [18], [191 [20] (see also [2T]). 

The factorization of the second possible term H ■ ii®S, namely the subleading 
jet functions, was proven recently in [22]. The subleading jet functions are zero at 
tree level and they were calculated explicitly at order as in [22J. As a result their 
contribution is suppressed by AgcD/^^b and a^. They are relevant for a high precision 
extraction of \Vuh\ using charmless and radiative inclusive B decays. Since they are 
part of the direct contribution, they reduce to local operators in the integrated rate. 

In summary, the three types of new perturbative calculations are at NNLO level 
in Us and/or Aqcd/'^6- They are almost at the theoretical limit, which implies that 
further improvement seems unlikely. Numerically the new perturbative corrections 
amount to about 1% correction for r(_B — )■ Xs7). 

3 Recent developments in ^ ^ Xg^' 
Non-Perturbative Aspects 

It was a common misconception in the field to assume that just like other inclusive 
B decays non-perturbative effects arise at 0(AQCj-)/m^) . While this is true for the 
Q77 — Qi-y contribution, it is no longer true when other operators are included. Over 
the years there were hints that "not all is well" in the study of the Q^g — Qsg [23], l2l] 
contribution and the Qi — Q^^ contribution [25] [26] [27] [28]. In particular in [26] it 
was stated, without a proof, that "There is no OPE that allows one to parametrize 
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non-perturbative effects from the photon couphng to hght quarks in terms of B meson 
matrix elements of local operators." Despite this statement, these effects were thought 
to be under control or small. There was never a systematic study of these effects. 
In fact, the uncertainty from Qj^ — Qsg p9] was largely misseqj. The conclusion 
of |29] was that non-perturbative corrections to the integrated rate arise already at 
^(AqcdM). 

A systematic study of these non-perturbative effects was performed recently in 
[5], which established the factorization formula ([1]) for the spectrum in the endpoint 
region. The new ingredient is the resolved photon contribution for which the photon 
does not couple directly to the weak vertex. In the integrated rate the resolved photon 
contribution does not reduce to matrix elements of local operators. Schematically, one 
finds terms of the form AF ~ J /i, where J is calculable in perturbation theory and 
h is a non-local matrix element. More specifically, for a photon energy cut > Eq 
define 

OPE 



OjlOPE 



where A = mi, — 2Eq and r(_E'o)|oPE is the "older" calculation. Assuming A ^ Aqcd 



-r (A\ '^i(^) Ai7(mg/mb,/i) Csgjfi) , . A^f '(^) 



rrih 



TTlb 



+ 



■ A88(A,^) Cpasifi) A A 

Anasifi) In — 

97r mi, ms 



+ 



(3) 



where model independently, 



Ai7 — ,/i 



Ar(/^) 

A88(A,/i) 



Cc Re 



Re 



oo ^1 
°° dUi 



12ml J 



mbUJi 

dU2 



/ll7(wi,/i) 



^ui + ie UJ2 - le 



uji + ie 



U2 - le 



/;,78^^^(wi,a;2,/i) , 

88 (A,c<;i,c<;2,yu) 



%Afln^ 
87r2 V A 



(4) 



F arises from a charm quark loop and its explicit form can be found in [5]. Auv 
is introduced to regularize the convolution integrals, but A88(A,/i) is independent 
of Auy. The exact definition of /ijj appears in f5], but schematically they are (F.T. 



^Spectator effects from Qy-y — Q^g were considered in [30] but tliey were estimated in a model 
dependent way that underestimated tfieir effects. 
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denotes the Fourier transform), 



/i88(wi,W2) F.T. of {B\b{0)---s{un)s{rn)---b{0)\B) 
hniioi) F.T. of {B\b{0)---G{sn)---b{0)\B) 

hjsiuj^,U2) F.T. of (5|6(0)---6(0)^ e,g(rn)---g(sn)|S). (5) 

Estimating the non-perturbative error for r(S — )■ Xg'j) reduces to estimating the size 
of Aij. Naively, for Aij ~ Aqcd ~ 0.5 GeV, the effects on the rate can be up to 30%. 
Fortunately, it is possible to constrain A17 and A^g'^"^. We now consider each of these 
parameters separately. 

3.1 Ai7 

The soft function hnlui) is an even function of ui and its normalization is 2A2 ~ 
0.24 GeV^. It might be natural to model hiY{uJi) as an exponential or a Gaussian 
which leads to an estimate of —10 MeV < A17 < MeV. But this is not a conservative 
bound, since hij{ui) does not have to be positive. Using the models of [5J, one finds 
—60 MeV < A17 < 25 MeV. A17 is smaller than the naive power counting estimate, 
since part of the contribution is already included in r(i?o)|oPE [23 ESI 123 ESj. It is 
also suppressed since F is peaked around 1 GeV, where hn{ui) is already becoming 
small. 

3.2 A^P"" 

In order to estimate Ayg'^'^ one can try and use one of the two methods. The first is to 
use the Vacuum Insertion Approximation(VIA). By using Fierz transformation one 
can write hjs{^i,^2) as a product of two B- meson light cone distribution amplitudes 
(LCDAs). A^^^'^ depends on the LCDA's inverse moment A^. One has A7g^'^|^j^ e 
Cspec [-386 MeV, -35 MeV], where Cspec = -1/3 for 5° or 5°, and e^pec = 2/3 for B^. 
The second approach is to assume SU (3) flavor symmetry. In this limit, Ayg'^^ is deter- 
mined by the isospin asymmetry Aq- in i? — > Xs'~f pi] . This asymmetry was measured 
by BaBar using two different methods. The naive average is Aq- = (— 1.3±5.9)%. In- 
cluding 30% SU{3) flavor breaking gives, Af^''^ ^ -4.5 GeV (cspec ± 0.05) Aq-. Notice 
that in the flavor averaged rate one has Cgpf^ = 1/6 which is effectively a suppression 
factor. 

3.3 Ass 

In this case one models Agg(A,/i) by Agg(A,/i) ^ A(/i) where A(yu) is positive and 
taken to be in the range < A(yLi) < 1 GeV. 
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3.4 Total Uncertainty 

Using the above values one has ^e\-^'j ^ [—1.7, +4.0] % and -^^Igg G [—0.3, +1.9] %. 

For J^e\-,s one has J^e\^1^ e [-2.8, -0.3] % or J-^l'^ e [-4.4, +5.6] % at (95% CL). 
"Scanning" over the ranges one has —4.8% < J-e(A) < +5.6% using VIA for A^^'^'^ 
or —6.4% < J^£;(A) < +11.5% using A^^'^'^ from Ao_. The last value reflects the 
large error on Aq-. But even if the error on Aq- was zero, one would still have 
-4.0% < J'e(A) < +4.8% in this "ideal case". 

In conclusion, one finds a total uncertainty of about 5%. This is also the un- 
certainty from a previous estimate based on Qj^ — Q^g alone [29], with an extra 
50% deviation from VIA. This value is also the uncertainty usually assigned to non- 
perturbative effects in various SM predictions. While numerically the prediction has 
not changed, it is now based on a much stronger theoretical basis. Still to be done are 
the analysis of the resolved photon contribution to the CP asymmetry, the spectrum, 
and the non-perturbative parameters extracted from it [32] . 

4 Comments on B Xd^y 

B — i- Xrf7 is analogous to 5 — )■ Xs7, where usually one only needs to replace V^^V^^ 
by VgbV*^. One important difference is that for this decay is not CKM suppressed 
and as a result the — Qj^ contribution is not CKM suppressed either . In [28] this 
contribution was estimated to scale as (9(Aqcd/''^^6), although it was not calculated 
explicitly. In [5] the contribution of — Qi-y to the CP averaged rate was calculated 

explicitly. It is given in terms of a convolution of P — and hi-j{(jji), i.e. the first 

term in the first line of (HI). Since the former is odd and the latter is even, the 
convolution vanishes. This removes the largest source uncertainty and makes B — )• 
Xff) as theoretically clean as S — i- Xg'^ [33] . 

5 Comments on B Xs 

For a more detailed discussion oi B ^ Xs see E. Lunghi talk at CKM 2008 [34J. 
Here we limit ourselves to a short review of recent developments. In the region of 
low g^, i.e. e [1...6] GeV^ and when one introduces a cut on the invariant mass of 
Xg, namely mx < m™*, dVi of B ^ Xg factorizes similarly to dTn of B Xg'y 
in the endpoint region, see (|2]). Here dTi, with i = T,A,L, corresponds to angular 
decomposition of the triple spectrum (called Hi in [35]). 

In [35j the contribution of the "primary" SSF to dTi was calculated. The primary 
SSF are defined to be the SSF that also contribute to 5 — )■ XJ v in the endpoint 
region. The authors of [53] find sizable corrections of the order 5% to 10% from the 
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primary SSF corrections. These cause a shift of ~ — O.OSGeV^ to — O.lGeV^ in the 
zero of the forward-backward asymmetry. 

In [36] the two-loop calculation of the hard functions was performed. The authors 
of [36] have observed a significant shift in the zero of the forward-backward asymmetry 
going from NLO to NNLO. Including the primary SSF contribution, they find the 
location of the zero to be at ql = (3.34 ... 3.40)^[J;| GeV^ for m^^* = (2.0 ... 1.8) GeV. 

Finally, following the completed analysis for r(i? — t- Xg'-f) one should ask what is 
the effect from "non-primary" SSF? For example, from soft gluon attachments to the 
charm-loop diagrams. This point was also stressed in [36] and requires further study. 

6 Summary and Outlook 

Inclusive Radiative B decays is a mature field. New factorization formula for photon 
spectrum in the endpoint region has been established in |5], which includes apart from 
the direct photon contribution familiar from B — ?■ Xu Iv, also new resolved photon 
contribution. The resolved photon contribution leads to (9(AQCD/'"^fe) corrections to 
T{B — )■ Xs 7). From a systematic study of these effects, one finds an irreducible error 
of ~ 5% on r(i? — >■ This non-perturbative error is the largest of the errors of 

the SM prediction, and therefore there is no prospect for reducing the total theoretical 
error below the 5% level. In the near future we can expect the perturbative error to 
be reduced below the non-perturbative and the experimental errors, and to have an 
analysis of the resolved photon contribution to the CP asymmetry in i? — )■ ^^7, the 
spectrum, and the non-perturbative parameters extracted from it . For B ^ 
the effects of the resolved photon contribution are yet to be calculated. 

Beyond that, further theoretical improvement seems unlikely. From the pertur- 
bative side, improvement beyond NNLO seems almost impossible and likely to be 
unjustified considering the non-perturbative and experimental errors. From the non- 
perturbative side, we are facing irreducible hadronic uncertainties. This implies that 
in the very near future we will have the definitive SM theoretical predictions for 
inclusive radiative B decays in years to come. 
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